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ABSTRACT 

Strong magnetic fields in relativistic stars can be a cause of crust frac- 
turing, resulting in the excitation of global torsional oscillations. Such os- 
cillations could become observable in gravitational waves or in high-energy 
radiation, thus becoming a tool for probing the equation of state of relativis- 
tic stars. As the eigenfrequency of torsional oscillation modes is affected by 
the presence of a strong magnetic field, we study torsional modes in mag- 
netized relativistic stars. We derive the linearized perturbation equations 
that govern torsional oscillations coupled to the oscillations of a magnetic 
field, when variations in the metric are neglected (Cowling approximation). 
The oscillations are described by a single two-dimensional wave equation, 
which can be solved as a boundary value problem to obtain eigenfrequcn- 
cies. We find that in the non-magnetized case, typical oscillation periods of 
the fundamental 1=2 torsional modes can be nearly a factor of two larger 
for relativistic stars than previously computed in the Newtonian limit. For 
magnetized stars, we show that the influence of the magnetic field is highly 
dependent on the assumed magnetic field configuration and simple estimates 
obtained previously in the literature cannot be used for identifying normal 
modes observationally. 
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1 INTRODUCTION 

Although the interior composition of relativistic stars is currently still very uncertain, the properties of their solid 
crust have been studied extensively (Ruderman 1968, see Pethick & Ravenhall 1999 for a recent review). It has 
been suggested that the presence of a strong magnetic field in a secularly evolving star can cause the crust to 
fracture, exciting shear waves and leading to the phenomenon of soft gamma repeaters (SGRs) (Cheng et al. 
1996, see Thompson 2000 for a recent review of SGRs). In such a scenario, the oscillation modes that appear 
most favored for excitation are the low-order torsional oscillations of the crust (Duncan 1998) . Torsional oscillations 
differ from spheroidal shear oscillations, in that they are predominantly divergence-free, toroidal velocity oscillations, 
accompanied by only small density oscillations in the presence of rotation or a magnetic field. Although torsional 
oscillations of neutron star crusts have a strong potential for observational detection, the study of these modes in 
the literature has been limited, so far, to a few representative models of the neutron star equilibrium and crust 
structure, while the magnetic field has only been taken into account in an approximate, Newtonian framework. 

Torsional modes may be favored for excitation (compared to spheroidal shear modes) during a starquake, 
because the restoring force is mainly due to the relatively weak Coulomb forces of the crustal ions (thus requiring 
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much less energy than compressional oscillations) and because the lowest-order torsional oscillation frequency has a 
relatively long period, which implies a slow damping rate (Duncan 1998). Although no oscillation modes in neutron 
stars have been detected to date, there is evidence for a 23 ms periodicity in the initial pulse of the 1979 March 5 
event (Barat et al. 1983). If confirmed, then torsional modes could become the first normal modes in neutron stars to 
be detected. As the frequencies of the torsional modes are significantly affected in the presence of a strong magnetic 
field, accurate frequencies for these modes for various magnetic field strengths and magnetic field configurations are 
needed, if SGRs are indeed strongly magnetized compact stars. 

In the Newtonian framework, torsional modes were first studied by Alterman, Jarosh & Pekeris (1959), applied 
to oscillations of the Earth's solid crust. That torsional modes could also be excited in neutron star crusts was first 
proposed by Ruderman (1968) and studied by Hansen & Cioffi (1980), Van Horn (1980), McDermott et al. 1985 and 
McDermott, Van Horn & Hansen (1988). Torsional modes of slowly rotating neutron stars have been studied by 
Lee & Strohmayer (1996). In general relativity, the theory of torsional oscillations was developed by Schumaker & 
Thorne (1983), applying the general-relativistic theory of elasticity by Carter & Quintana (1972) (for a discussion, 
see Priou 1992 and references therein). The only numerical computation of torsional modes for realistic neutron 
star crusts in general relativity has been presented by Leins (1994). 

In the present work, we initiate the study of torsional oscillations in relativistic stars possessing a strong 
magnetic field. Our theoretical description is based on the study of wave propagation in hydromagnetic media 
in general relativity by Papadopoulos & Esposito (1982). Previously, the magnetic field effect on the torsional 
modes has been studied in the Newtonian limit by Carroll et al. (1986) and by Nasiri & Sobuti (1989) (see also 
Duncan, 1998). Here, we derive the linear perturbation equations governing torsional oscillations in a magnetized 
star in the relativistic Cowling approximation (Cowling 1941; McDermott, Van Horn & Scholl, 1983), neglecting 
the deformation of the equilibrium structure, due to the presence of a magnetic field. For a general axisymmetric 
magnetic field configuration, the perturbation equations are two-dimensional. Simplified, one-dimensional equations 
are derived for a special case of the magnetic field configuration. In future work, we plan to include the deformation 
of the neutron star structure and study more general magnetic field configurations by solving the two-dimensional 
eigenvalue problem for various magnetic field configurations. 



2 THE EQUILIBRIUM CONFIGURATION 

The equilibrium structure of a magnetized relativistic star is non-spherical, due to the deformation induced by 
magnetic field stresses. Self-consistent models of relativistic stars with a strong magnetic field have been constructed 
by Bocquet et al. (1995) (in the case of rapidly rotating stars) and Gupta et al. (1998) (in the slow-rotation limit). 
The spherical structure of a nonrotating relativistic star is distorted significantly, only when the magnetic field 
becomes very large, exceeding ~ 10 14 G (for a magnetic dipole). Although this distortion will certainly modify 
the eigenfunctions of various oscillation modes, we expect this effect to be smaller than the direct influence of 
the magnetic field on the eigenfunctions and eigenfrequencies of oscillations, unless the magnetic field becomes 
exceedingly large. Thus, here we will neglect the distortion induced by the magnetic field and will assume the 
equilibrium structure to be that of a spherical relativistic star, the space-time of which is described by the metric 

ds 2 = - e 2 *M dt 2 + e 2A ^dr 2 + r 2 (d6 2 + sin 2 6dtf) , (1) 

where $(r) and A(r) are functions of the radial coordinate only. The equilibrium structure is obtained in full general 
relativity (see e.g. Schumaker & Thorne 1983), assuming a perfect-fluid equation of state. Although our purpose is 
to study oscillations in the solid crust of the star, the bulk properties and local pressure and density profiles in the 
crust can be obtained accurately, by assuming a perfect fluid. We assume, further, that the crust has an isotropic 
shear modulus. The shear modulus only appears in the perturbed configuration. 

In order to derive the perturbation equations for a magnetized relativistic star, the equilibrium stress-energy 
tensor is assumed to consist of two parts: 

rpab rpab(pi) _j_ rpab(M) 

where 

T aft(pf) = (e + p)u a u b +pg ab , (3) 



© RAS, MNRAS 000, 



Torsional Oscillations of Magnetized Relativistic Stars 3 



and 

T ab(M) =R 2 u a u b + H^_ g ab _ R a H b (4) 

(see Papadopoulos & Esposito, 1982). Above, e and p are the energy density and pressure of the fluid and u a is the 
four-velocity of fluid elements 

* M ,0,0,0l , (5) 



u = e 



(with u a u a = —1), while H a are the components of the magnetic field (with H 2 — g ab H a H b ). The solid crust is 
assumed to be isotropic and perfectly elastic, so that in equilibrium there is no shear stress contribution in (^|). 
Here, and throughout the paper, we are using the same units as in Papadopoulos & Esposito (1982), i.e. we set 
c = G = 1 and normalize the magnetic field by multiplying H a by \^4jt, unless otherwise noted. 

The equations of motion of the fluid are obtained by the projection of the conservation of the stress-energy 
tensor onto the hypersurface normal to u a 

h c a T ab , b = 0, (6) 

(where h c a = g c a + u c u a and ";" denotes the covariant derivative compatible with the equilibrium metric g a b), 
which reads 

(e + p + H 2 )u b u a , b = - (p+^Pj h ab + h a c (H c H b ) ]b . (7) 

For perfectly conducting hydromagnetic media, Maxwell's equations F ab -b = (where F ab is the electromagnetic 
tensor) take the form 

(u a H b -u b H a ), a = 0, (8) 
which leads to the following form of the magnetic induction equation 

H a . b u b = (a a b + tu a b - j^8 a b@^ H b + H b Ub- c u c u a , (9) 

where = u b ; b is the expansion of the fluid, ui a b is the twist tensor and a a b is the shear tensor. We assume a zero 
net electrical surface charge, so that there is no equilibrium electric field. The assumption of infinite conductivity 
is justified as the magnetic field diffusion timescale is several orders of magnitude larger than the typical timescale 
for torsional oscillations. 



3 THE PERTURBED CONFIGURATION 

The various oscillation modes that can be present in a star are studied by linearizing the equations governing the 
equilibrium configuration and assuming a harmonic time-dependence. Here, we are only interested in the torsional 
modes of the crust, which do not generate significant variations in the gravitational field of the star. In addition, 
the low-order quadrupole torsional mode has a long period of the order of 20-40ms. Thus, such an oscillation can 
be described well by the slow-motion approximation (Schumaker & Thorne 1983). It follows that one can describe 
the low-order torsional oscillations of the crust by neglecting the variations in the metric, i.e. by setting Sg a b = 0, 
which is usually called the relativistic Cowling approximation (McDermott et al. 1983) (this approximation was 
also used in the description of hydromagnetic waves in general relativity by Papadopoulos & Esposito, 1982). This 
approximation yields the real part of the eigenfrequencies with a typical accuracy of a few percent, as can be 
deduced from the results in Leins (1994) and we will employ it throughout this paper. Once the eigenfunction and 
the real part of the eigenfrequency are computed, estimates for the damping rate of torsional oscillations due to 
gravitational and electromagnetic wave emission and due to viscous dissipation can be obtained in a perturbative 
way (see e.g. McDermott et al. 1988) , since the imaginary part of the eigenfrequency is much smaller than the real 
part. 

The linearized version of the equations of motion (0) is 
(e + p + H 2 ^ 5u a -bU = — (Se + Sp + 2H C 8H C ) u a - b u b + (u a 8u c + 5u a u c ) 
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-(e + p + H 2 ) u a , b Su b + h\ [H c SH b + SH c H b - g cb (Sp + H d 8H d )] ;f> 
-h a c 8T ch f. (10) 
In ([H]), the perturbation in the shear stress tensor is 

STjV = ~2^8S ab , (11) 

where SS a b is the perturbation in the strain tensor and fi is the isotropic shear modulus. The linearized version of 
the magnetic induction equation (^) is 

{SH a )-, c u c = -H a ; C 6u c + h ac H d {6u c ); d + h ac u C ; d 6H d + (u a Su c + u c Su a )H d u c . d - 65H a - H a S6 

—u a H° [(8u c )- l bU b + u c[ bSu b ~^ — (H c Su a + u a SH c ) u c - b u b . (12) 

Equations @ and 0) are the complete set of equations that govern the perturbations of magnetized relativistic 
stars in the Cowling approximation. 



4 THE EIGENVALUE PROBLEM 

When the distortion of the equilibrium structure, due to the magnetic field, is ignored, the spherical symmetry of the 
unperturbed star allows for the oscillations to be decoupled into modes of definite spherical-harmonic indices (I, m) 
and definite parity. Here, we investigate pure torsional oscillations, which are the normal modes of odd (magnetic- 
type) parity, tv = (— 1) !+1 (Regge & Wheeler, 1957). In spherical symmetry, modes with fixed I but different m yield 
the same frequency, thus, we will specialize to the case of m = only. This means that, we will not study (at this 
point) the mode-splitting caused by the magnetic field. 

When expanded in vector spherical harmonics of definite I and m — 0, the odd-parity perturbation in the 
four-velocity can be written as (see Schumaker & Thorne, 1983) 



Su = 
<5i/ = 
where 



dt 



1 dPi(cost 



(13) 
(14) 

(15) 



sin 6 d6 

and Y(r,t) is the angular displacement of the oscillating stellar material and Pi(cos9) is the Legendre polynomial 
of order I. For odd-parity perturbations in a spherical background: 8u = Su r = 5e = 0. 

We assume a harmonic time-dependence of all perturbed variables, as in Y(r,t) = y(r)e 1 "*, where u is the 
mode frequency (from now on we drop the time dependence from all perturbed variables). The </>— component of 
the perturbed equations of motion ([To]) becomes 



iu (e + p + H )e 5u 



5H r 
+SH 



H^, r + H 



$.r + A.r + ~ 

r 







+ 5H 9 


1 r 2 sin 2 6 





#% + 3cot OH* 



2 cotOH 6 + -H r 



H r {A, r + <f,r 

+H 9 [5H\ r + 5H\ e + 5H 4, + H B 5H" 

( -!■ • •!> , • \.) ^/ -3c.it OST 



+ H r 8H' t ' 
+ H 4 '8H' 1 ' 



,4> ' 



-ST 



MS) 



ST 



8<p(S) 



r 2 sin 9 

6I0(S) 



H 9 5Hg 



where 

n(S) 
L r<j> 

and 



_ $ BY ^ 
dr ' 



(16) 



(17) 



ST^' = -nr sin 



Oe^Y- 



30 



(18) 
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are the only non-vanishing components of the perturbed shear stress tensor. The oscillating magnetic field compo- 
nents are given algebraically as 



5H r 
5H e 



SH 9 = ^- (H r 

iu> 



, r + H B 8u' 1 ' 



(19) 
(20) 

(21) 



Equation ( |l6| ) represents a two-dimensional boundary value problem for the eigenfrequency lj and the eigenfunction 
Y(r). If we substitute (|l7|)-(^l|) in (|l^), we arrive at our final expression 



(| + 2)(i-l)/i^ r Y- 



h 5>,r — A >r ) fj, + H tJ . 



2(*-A) 



| fo* H r (W 9 cot + ~iT + H r (A, r + $ r )) {-$, r Y + Y r ) + H r IT r (-$, r Y + Y, r ) 



+i? e iy ; 9 (-$, r y + f,.) + (fry (-$, r Y + y,r 



1 1® 



2iTY H cot 



+H e H%Y 



+ b% e (H e ) 2 Y 




H r (H e r Y + H e Y r ) + H°H r (-$, T Y + Y r ) 



(22) 



In the absence of a magnetic field, separation of variables allows the above equation to be reduced to a one- 
dimensional equation with respect to the radial coordinate. This is not possible for a general magnetic field con- 
figuration. However, as we will show in Section ^, there exist special cases of the magnetic field, for which the 
problem becomes one-dimensional. As the numerical solution of the two-dimensional problem is not a trivial task, 
such one- dimensional cases are very useful for obtaining initial order-of-magnitude estimates for the influence of a 
strong magnetic field on the torsional oscillations. 



5 ANALYTIC ESTIMATES 

If one assumes a uniform density star with uniform shear modulus fi in the Newtonian limit, one can derive a simple 
analytic estimate for the influence of the magnetic field on the period of torsional oscillations. In the nonmagnetized 
case, the period of torsional oscillations is given analytically as 

P=^, (23) 

where v s = \ ffjjp is the speed of sound, p is the density of the crust, R is the star's radius and x n is a constant 
(see Schumaker and Thorne 1983). In the presence of a magnetic field B, the density p is replaced by p + B 2 /4ir. If 
one also assumes that the shear modulus /i is augmented by the magnetic field tension B 2 /Ait, one easily obtains 



p ^ i + ( B /lr (24) 

where va = B 2 /4np is the Alfven speed, B M = {A-np) 1 ^ 2 and Pq is the oscillation period for a nonmagnetized star. 
This can be rewritten as 



P = Po\l V,T D 7,r?2 - + (B/^) 2 ]" 1/2 , (25) 



1 + vUB/B^ 

1 + [b/b,y 

when v s « 1 (which was also considered by Duncan, 1998). The above estimate implies that the main influence of 
the magnetic field on the oscillation period is through the magnetic field tension. Based on a power-law fit to the 
deep-crust equilibrium composition by Negele and Vautherin (1973), Duncan (1998) estimates _B M = 4 x 10 15 p ( 1 4 4 G, 
where pi4 = p/10 14 grcm -3 . This then yields a significant decrease in the oscillation period at a few times 10 15 G. 
As we will show in the next Section, this simple estimate is misleading, as it completely ignores the magnetic field 
configuration. In practice, the restoring force for the torsional oscillations will not be augmented by the magnetic 
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field tension equally throughout the star, but the influence of the magnetic field on the oscillation will be highly 
dependent on the correlation between the multipoles of the magnetic field configuration and the quadrupole toroidal 
velocity field of the torsional mode. For example, as we are showing in the next section, the radial component of the 
magnetic field has a much weaker influence on the torsional modes than other components. Thus, simple estimates as 
the one made above, are not sufficiently accurate for obtaining quantitative results that could be used for identifying 
normal modes in observations. 



6 A SPECIAL CASE FOR THE MAGNETIC FIELD 



In order to arrive at a one-dimensional boundary value problem, we consider the following magnetic field configu- 
ration as a toy-model: 



H 9 = 0, 

For this case, Maxwell's equations reduce to 
2^ 



H 



A, 



0. 



The solutions admitted by (j2S 

e~ A 
H = — 5-/*o, 



are of the form 



(26) 
(27) 

(28) 
(29) 



where /io is a constant. Although this form of the magnetic field cannot be considered realistic, we will show that 
this toy-model can be used for obtaining first estimates of the influence of the magnetic field on the oscillation 
frequencies. In our derivation, we assume that the magnetic field has the above form only inside the solid crust and 
not throughout the star. 

Using (^), one can reduce the two-dimensional boundary value problem (H) to the following one-dimensional 
form: 



(e+p + H 2 ) 



4 A- 

r e 



•y 



{[,i+( e A /r) 2 ] e *-Vr r } +(A, r 

- {(I + 2)0 - l)p + r 2 [($ r ) 2 + $,„.) 



i \ 4>+A 4 

<P.rje r 



(HT Y r 



(HI 2 } 



S+A 2 V 

e r Y. 



(30) 



This equation is a generalization of equation (66a) in Schumaker & Thorne (1983) to the case of the magnetic field 
configuration assumed above, in the Cowling approximation. We note that both the influence of the equilibrium 
magnetic field as well as the coupling between torsional and magnetic field oscillations are included in our derivation. 

Comparing eqn. (BO) to the non-magnetized case, on sees that the magnetic field enters through different terms. 
First, on the l.h.s. of (3C), the total energy density of the hydromagnetic medium becomes e + p + H 2 , as one would 
expect. On the r.h.s., an important contribution of the magnetic field is the augmentation of the shear modulus /i (in 
the term involving the second spatial derivative of Y) by the magnetic field term (e A H r ) 2 . This is the only magnetic 
field term on the r.h.s. that survives in the Newtonian limit, which is in agreement with the Newtonian results in 
Carroll et al. (1986) (see also Duncan, 1998). The other two contributions of the magnetic field on the r.h.s. are 
purely relativistic effects. Carroll et al. (1986) only looked at a cylindrical model of the neutron star crust near the 
magnetic polar cap, with a uniform magnetic field normal to the surface of the crust. The one-dimensional boundary 
value problem considered in this section is a generalization of Carroll et al.'s cylindrical Newtonian treatment to a 
spherical geometry and to general relativity. 

Eigenvalues of the torsional modes in magnetized stars can be obtained by solving ( po| ) as an eigenvalue problem 
with zero-traction boundary conditions at the base of the crust and at the surface of the star. For the oscillations 
of the magnetic field, we assume approximate boundary conditions, neglecting the oscillations in the exterior of the 
solid crust. For the numerical solution of the above second-order eigenvalue problem, we introduced new variables 
Yi and Y-z (see Carroll et al. (1986) and Leins (1994) for similar definitions): 



Yi 



Y, 



Y 



4tt 



dr 



(31) 
(32) 
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Figure 1. Oscillation period of the fundamental 1=2 torsional mode as a function of the magnetic field strength, for a 1.4Mq 
relativistic star and for the simple magnetic field configuration considered in Section tol The influence of the magnetic field 
is opposite to that predicted by simple estimates (see Section ph, which demonstrates that the oscillation period is highly 
dependent on the given magnetic field configuration. 



(note that we have restored the factor of 4/7T dividing the magnetic field). This has the advantages of avoiding the 
numerical evaluation of the derivative of the shear modulus \i (which is known in a tabular form only, for realistic 
equations of state) and that Y 2 vanishes at the boundaries when zero-traction is assumed. With this definition, the 
problem of finding the eigenfrequencies uj of torsional modes reduces to solving the system of equations 



dYi 
dr 

dY 2 
dr 



Y 2 



+ 



477 

4 (A 
r 



(33) 



$, r )(e A H r 



4tth + (e A H r ) 2 
(l + 2)(l-l)fi + r 2 



Y 2 

( — 

\~dr 



2 d 2 ® 
dr 2 



(H r ) 



4tt 



P + 



EL 

4tt 



v l5 



(34) 



with boundary conditions Y 2 — at the surface and the base of the crust and Y\ being any finite number at the 
base of the crust. 

We construct an equilibrium TOV model, using the Wiringa et al. (1988) equation of state (WFF3) matched to 
the Negele and Vautherin (1973) equation of state at low densities. The star has a gravitational mass of M=1.4 Mq 
and a radius of R = 10.84 km. We assume that the crust extends from energy densities of 2.4 x 10 14 gr/cm~ 3 down 
to energy densities of 5 x 10 10 gr/cm -3 , which corresponds to radii from 9.86km to 10.65km (we truncate the crust 
at very low densities to avoid numerical difficulties - this does not significantly affect the obtained eigenfrequencies) . 
We note that, apart from the work of Leins (1994), the frequencies of even the non-magnetized case have not been 
computed for realistic equations of state in general relativity. In all previous Newtonian studies, the eigenfrequency 
of the fundamental quadrupole torsional mode has been found to be close to 20ms. In our example, this frequency 
is 37ms. This shows that the structure of a general relativistic star and the relativistic perturbation equations yield 
torsional frequencies very different from those predicted by the Newtonian equations. For example, part of the 
difference between the Newtonian and relativistic eigenfrequencies is due to the relativistic term e 2 *, multiplying fi 
in equation (^). For our model, e 2 * = 0.59 at the base of the crust. This effect increases the fundamental torsional 
period by a factor of roughly 1/V0.59 ~ 1.3. The remaining difference between our obtained 37ms fundamental 
period and the 17.32 ms fundamental period for the most massive model considered in McDermott et al. (1988) 
(model NS13T8) can be explained as follows: the model considered in McDermott et al. (1988) has a mass of 1.326 
Mq, similar to the mass of our present model, but the radius is only 7.853 km, (compared to a radius of 10.65km for 
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our model) . The period of the fundamental torsional mode has been shown to be roughly proportional to the star's 
radius (when the crust is thin) in Hansen and Cioffi (1980), who estimate the I — 2 period as 24.5(i?/10km)ms. 
Using the radius of our model and correcting for general relativity, the above estimate becomes 34.5ms, which is 
close to our numerically obtained result of 37ms. Thus, a general-relativistic version of Hansen & Cioffi's formula 
for the period of the fundamental torsional mode is 



which should hold for 1.4M© models constructed with various realistic equations of state, with the possible exception 
of strange star models with very thin crusts. 

Figure 1 shows the period (in ms) of the fundamental quadrupole torsional mode as a function of the magnetic 
field B = H r . Due to the relativistic factor which counteracts the shear modulus in (Q), the period is increased by 
the magnetic field for values B > 2 x 10 15 G. The change in the period does not follow the dependence predicted by 
the simple estimate (|if]) , which predicted that the magnetic field decreases the oscillation period at a few times 
10 15 G. We conclude that the eigenfrequency of torsional modes in the presence of a strong magnetic field cannot 
be simply estimated, assuming that the magnetic tension augments the shear modulus in restoring the oscillation, 
isotropically. Instead, the degree to which the magnetic field modifies the eigenfrequency is highly dependent on the 
magnetic field configuration and on the correlation of the latter with the toroidal velocity field of the oscillation. 

In our example, the magnetic field has only a radial component. A careful inspection of ( |34| ) and (^i|) reveals that 
the main term that determines the eigenfrequency, (I +2)(i — is not augmented by the radial component of the 
magnetic tension in the Newtonian limit, but only by the azimuthal component. Thus, for a general configuration 
of the magnetic field, having both radial and azimuthal components, the numerical results of Figure 1 can be 
considered as one limiting case of the influence of the magnetic field on the torsional modes, while the isotropic 
simple estimate ( p5| ) would be another limiting case. We expect that the frequencies of torsional modes for realistic 
magnetic field configurations will have values within these limits. 



7 DISCUSSION 

We initiate a study of torsional oscillations in relativistic stars possessing a strong magnetic field, based on the 
treatment of wave- propagation in hydromagnetic media in general relativity by Papadopoulos & Esposito (1982). 
We derive the linear perturbation equations governing torsional oscillations in a magnetized star in the relativistic 
Cowling approximation, neglecting the deformation of the equilibrium structure due to the presence of a magnetic 
field and variations of the metric. The perturbation equations are two-dimensional in the case of a general axisym- 
metric magnetic field configuration. Simplified, one-dimensional equations are derived for a limiting case of the 
magnetic field configuration. This allows first estimates for the change in the mode-frequencies, due to the magnetic 
field, to be obtained. The origin of this change is both due to the equilibrium energy density of the magnetic field 
and due to the back-reaction of the magnetic field oscillations on the oscillations of the solid crust. Our results 
are substantially different from a simple isotropic estimate, showing that the influence of the magnetic field on the 
frequency of torsional oscillations is highly dependent both on relativistic effects and the structure of the magnetic 



The study of torsional oscillations is motivated by the prospect that they could become observable in gravi- 
tational waves or in high-energy radiation after a crust fracture is initiated by a strong magnetic field. Work by 
Duncan (1998) suggests that in this scenario the torsional modes of the crust will be the dominant mode of oscil- 
lation, as they do not couple strongly to density variations. In the latter reference it is suggested that an observed 
periodicity of 21ms in the March 5 event could in fact be the signature of the fundamental quadrupole torsional 
oscillation of the neutron star's crust. For this to happen, the star must have a very strong magnetic field, since, as 
we show, this mode has a longer period of 37ms in a typical relativistic nonmagnetized star (contrary to the 20ms 
estimate derived in previous Newtonian studies). 

The correct identification of such periodicities with a specific normal mode requires the computation of mode- 
frequencies in the presence of a strong magnetic field in full general relativity. A successful identification will 
significantly constrain the properties of the high-energy equation of state in relativistic stars. As no other mode of 
oscillation has been identified in neutron stars, to date, torsional modes could be the oscillations to initiate the field 




(35) 



field. 
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of observational neutron-star asteroseismology. Especially a combination of mode-identifications obtained in high- 
energy radiation observations with mode-identifications in gravitational radiation detections (see Andersson and 
Kokkotas 1996, Kokkotas, Apostolatos and Andersson 2001) should yield invaluable information on the properties 
and internal composition of relativistic stars. 

In future work, we plan to present detailed numerical results for various realistic neutron star equations of 
state and for various axisymmetric magnetic field configurations, by numerically solving the full two-dimensional 
eigenvalue problem. 
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